The Fermi-Ulam Accelerator Model Under Scaling Analysis 
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The chaotic low energy region of the Fermi-Ulam simplified accelerator model is characterised by 
use of scaling analysis. It is shown that the average velocity and the roughness (variance of the 
average velocity) obey scaling functions with the same characteristic exponents. The formalism is 
widely applicable, including to billiards and to other chaotic systems. 
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Enrico Fermi [l| attempted to describe cosmic ray ac- 
celeration through a mechanism in which a charged parti- 
cle can be accelerated by collision with a time dependent 
magnetic field. His original model was later modified 
and studied in different versions, based on different ap- 
proaches, one of which is the well-known problem of a 
bouncing ball or Fermi-Ulam model (FUM) HH]. The 
main results for this version of the problem, consider- 
ing periodic oscillation, can be summarised as: (i) it is 
described via the formalism of an area-preserving map; 
(ii) it presents a set of invariant spanning curves in the 
phase space for high energy; (iii) a set of KAM islands 
surrounded by a chaotic sea can be observed in the low 
energy regime; and (iv) small chaotic regions limited by 
two different invariant spanning curves can be observed 
at intermediate energies. A related version of this prob- 
lem in a gravitational field, sometimes referred to as a 
bouncer 4], presents a property, in contradistinction to 
the bouncing ball problem, that depending on both ini- 
tial conditions and control parameters, the particle has 
unlimited energy gain i.e. the basic condition needed for 
Fermi acceleration. The difference between these appar- 
ently very similar models was latter clarified by Lichten- 
berg et al |5J. The quantum problems corresponding to 
both the FUM and bouncer models have also been inves- 
tigated la, 01 Q ■ The special interest in studying these 
one-dimensional classical systems arises because they al- 
low direct comparison of theoretical predictions with ex- 
perimental results 9,10]. Even more, the formalism used 
in the characterisation of such models can immediately 
be extended to the so-called billiards class of problems 

[HIIllil. 

In this Letter, we characterise the average velocity, and 
its variance which we will refer to as roughness^ within 
the chaotic sea of the phase space using scaling functions. 
One of our tools, the roughness, is an extension of the 
formalism used to characterise rough surfaces [iJI which, 
as we will show, is immediately applicable to chaotic or- 
bits in the problem of time dependent potential wells 
111 111 Elll as well as to billiards problems HI E^. 
This scaling scenario represents the first characterisation 
of the integrability-chaos transition in this problem and 
should be applicable to several billiard problems. The 
formalism may therefore prove useful in characterising 



classes of universality. 

Let us describe the system and how to characterise its 
dynamical evolution. It consists of a classical particle 
bouncing between two rigid walls, one of which is fixed; 
the other moves periodically in time with a normalised 
amplitude e. We will describe the system using a map 
T{Vn,4>n) = (Ki-i-i,^n-i-i) which givcs the new velocity 
of the particle and the corresponding phase of the mov- 
ing wall immediately after the particle suffers a collision 
with it. We will use a simplification 21] in our descrip- 
tion of this problem: we will suppose that both walls 
are fixed but that, when the particle suffers a collision 
with one of the walls, it exchanges momentum as if the 
wall were moving. This simplification carries the huge 
advantage of allowing us substantially to speed up our 
numerical simulations compared with the full model. It 
is usefully applicable because the main dynamical prop- 
erties of the system are preserved under such conditions. 
Incorporating this simplification in the model and using 
dimensionless variables, the map is written as Q 
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0n + T7- mod27r 



(1) 



The term 2/y„ specifies the length of time during 
which the particle travels between collisions, while 
— 2esin((/)„+i) gives the corresponding fraction of veloc- 
ity gained or lost in the collision. The modulus function 
is introduced to avoid the particle leaving the region be- 
tween the walls. We stress that the approximation of 
using the simplified FUM is valid in the limit of small 
e. So, the transition from integrability (e = 0) to chaos 
(e y^ 0), characterising the birth of the chaotic sea, can 
be well described. 

We will concentrate on the scaling behaviour present 
in the chaotic sea. We investigate the evolution of the 
velocity averaged in M initial phases, namely 



1 *^ 



(2) 



where Vq is the initial velocity and j refers to a sarnple 
of the ensemble. In order to define the roughness |14l |. 
we first consider the average of velocity over the orbit 



generated from one initial phase 
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Wc then evaluate the interface width around this aver- 
aged velocity. Finally, the roughness is defined by con- 
sidering an ensemble of M diff'erent initial phases: 



w(n,e,Vb) 



1 *^ 

-y 

M '^ 



TT^. , 



V^j{n,e,VG)-VAn,e,V(,) 



Fig-fflshows the behaviour of the roughness for two diff'er- 
ent control parameters. We can see in Fig. Q^a) that the 
roughness grows for small iteration number n and then 
saturates at large n. The change from growth to satura- 
tion is characterised by a crossover iteration number rix. 
We can also see that different values of the control pa- 
rameter e generate different behaviours for short n. This 
indicates that n is not a scale variable. However, it turns 
out that the transformation n -^ ne^ coalesces the curves 
for small n as illustrated in Fig.^b). We therefore infer: 
(i) that for n ^ n^ the roughness grows according to 

a;(ne2,e,Vo)(x(ne2)/3 , (4) 
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FIG. 1: (a) Behaviour of the roughness lj as a function of the 
iteration number n. (b) Behaviour of lj as a function of ne^. 
Both curves were derived from an ensemble average of 50, 000 
different initial conditions starting with Vb ~ 0. 
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FIG. 2: (a) Plot of LJsat against the control parameter e. (b) 
The crossover iteration number n^ as a function of e. 



where exponent /? is called the growth exponent; (ii) that, 
as the iteration number increases, for n ^ n^., the rough- 
ness reaches a saturation regime that is describable as 



a;sat(e) oc e° 



(5) 



where a is the roughening exponent; and (iii) that the 
crossover iteration number rix marking the approach to 
saturation is 



nx{e,Vo) oc e^ 



(6) 



where z is called the dynamical exponent. With these 
initial suppositions, we can now describe the roughness 
formally in terms of a scaling function. 



iu{ne^, e, Vq) = luj{l''ne'' ,Fe, FVq) , 



(7) 



where I is the scaling factor and a, 6, and c are referred 
to as scaling dimensions. It is important to stress that 
these scaling dimensions o, b and c must be related to 
the characteristic exponents a, (3 and z. All of the above 
discussion is also valid for the average velocity V . To 
relate the exponent (3 to the scaling dimensions, we chose 
/ = (ne^)^~. This allows us to rewrite iQ as 



The function oji ~ cj I 1, (ne^) <»£, (ne^) « Vq 1 is sup- 
posed constant for n <^ Ux- Comparing equations © 
and (01), we can however conclude that —-=13. Choos- 
ing Z = e^6 ^ we have that 



w(ne^,e, Vo) = e '•uj2{{n€^)€ 



■Vo) 



(9) 



where uj2 = (^ ((ne^)e~^, 1, e~^Vo) is assumed constant 
for n^ rix- Comparison of equations lO and shows 
that — -J = a. It is less straightforward to obtain the ex- 
ponent c. To do so, we use a result from a recent paper 
where two of us [23 utilised a connection with the well 
known Standard Model (SM) ^ to describe the position 
of the first invariant spanning curve (FISC) above the 
chaotic sea in the FUM. It was shown that the control 
parameter e could be related to a typical mean velocity on 
the FISC in the FUM to give an effective control param- 
eter Kcs = 4e/y*^ w 0.97 . . . which is gratifyingly close 
to the value of the control parameter Kc = 0.9716 ... at 
which the SM exhibits a transition from local to globally 
stochastic behaviour [23 ■ We can thus rewrite the effec- 
tive control parameter KcS in terms of scaled variables 
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FIG. 3: (a) Roughness evolution for different values of the 
control parameter e. (b) Collapse of the curves from (a) onto a 
universal curve. Both (a) and (b) were obtained using Vi) ~ 0. 
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FIG. 4: Behaviour of the average velocity V as a function of 
n for different values of e and Vq. (a) The original time series, 
(b) Collapse of the data onto a universal curve. 



as 



K, 



4(Z^ 



[I'VoY 



4e F 



(10) 



which implies that 6 — 2c = 0. Using our result for the 
exponent b, we find c = — ^. Note that all scaling di- 
mensions are therefore determined if we can obtain the 
exponents a and (3 numerically. The exponent a is ob- 
tained in the asymptotic limit of large iteration num- 
ber and it is independent of Vq . Fig. |2Ia) illustrates an 
attempt to characterise this exponent using the extrap- 
olated saturation roughness. Extrapolation is required 
because, even after lO'^n^; iterations, the roughness has 
still not quite reached saturation. From a power law fit, 
we obtain a = 0.512(3) « 1/2. We can thus rewrite Eq. 
as 

u;(ne2, e, Vo) = {ne^g [{ne^r^^e, {ne^r^V^] . (11) 

In order to obtain /?, we observe that we have two "time" 
scales in equation (|11|) . namely n'^ and n" and that the 
second one (n") is basically zero if we chose Vq « 0. 



Then we determine f3 from the short "time" behaviour 
{n ^ n'^). After averaging over different values of the 
control parameter e in the range e G [10^**, 10~^], we 
then obtain /3 = 0.496(6) « 1/2. Therefore, the scaling 
dimensions describing the scaling of the chaotic sea in 
the limit of small e are a ^ b = —2 and c = — 1. From 
Eqs. and JS)) we find that the scaling relation for the 
exponent z is z = a/ (3 — 2. Considering the previous 
values of both a and /3, we obtain that z = — 1. The 
exponent z can be also obtained numerically. Fig. I^b) 
shows the behaviour of the crossover iteration number n^. 
as function of the control parameter e. The power law 
fit gives us that z — —1.01(2), in good accord with the 
scaling result. The scaling for Vb ~ is demonstrated in 
Fig. 121 where the three different curves for the roughness 
in (a) are very well collapsed onto the universal curve seen 
in (b) when we normalise the quantities with a = b = —2. 
The case with initial velocity Vq y^ is better illus- 
trated by the average velocity (see Fig. 2J). Now, we 
must consider two "time" scales, namely n'^ oc 1/e and 
n" oc lo^/e^. From equation ^U^ (see also Ref. [23), 
the maximum initial velocity inside the chaotic sea is 
Vo,max ~ 2e^/^ implying that the second time scale has a 
maximum value of (n" '^ 4n^). So we observe that two 
different kinds of behaviour may occur, for n" < n'^ or 
"" ~ "a;- When Vb = 10~^, we have n^' ~ and we can 
see in Fig.^Ja) that the curves for e = 10"'' and e — 10"'^ 
show only two regimes: (1) a growth in power law for 
n <ti n'^ and (2) the saturation regime for n^ n'^. Con- 
sidering Vb = 10~^ and e = 10"'* we have that n" < n'^ 
and we can see for such curve in Fig.^fa) three regimes. 
For n ^ n", the average velocity is basically constant. 



When n" < n < n^ , the curve growth and begin to follow 
the curve of Vb = 10^^ and same e. In this window of 
n, we have a growth with a smaller effective exponent f3. 
Finally, for n '^ n'^ we have the saturation regime. It is 
shown in Fig 0Jb) that the collapse of the curves holds 
even for Vb ^ 0, implying that the inferred scaling form 
V^(ne^,e, Vb) with exponents a = b = —2 and c = —1 is 
also correct. 

In summary, we have characterised the average veloc- 
ity and its variance (roughness) in the chaotic sea in the 
simplified FUM by use of a scaling function. We show 
that the critical exponents /3, a and z are connected by 
a scaling relation. We emphasise that this behaviour 
is valid for small values of e and it is immediately ex- 
tendable to other average quantities. We have charac- 
terised, for the first time, the scaling appearing in the 
integrability^chaos transition (from e = to e ^ 0) of 
the FUM. The scaling scenario should also hold for bil- 
liard systems, so that this kind of formalism should be 
useful for characterising asymptotic properties in such 
problems. It should be possible to extend it to encompass 
time-dependent Hamiltonian systems and a huge class of 
other problems exhibiting chaotic behaviour. 

E. D. Leonel was supported by Conselho Nacional de 
Desenvolvimento Cientifico CNPq, from Brazil. The nu- 
merical results were obtained in the Centre for High 
Performance Computing in Lancaster University. The 
work was supported in part by the Engineering and 
Physical Sciences Research Council (UK). J. K. L. da 
Silva thanks to Conselho Nacional de Desenvolvimento 
Cientifico (CNPq) and Fundagao de Amparo a Pesquisa 
de Minas Gerais (Fapemig), Brazilian agencies. 



[1] E. Fermi, Phys. Rev. 75, 1169 (1949). [13] 

[2] A. J. Lichtenberg, M.A. Lieberman. Regular and Chaotic 

Dynamics{Appl. Matli. Sci. 38, Springer Verlag, New [14] 

York, 1992). 
[3] M.A. Lieberman and A. J. Lichtenberg. Phys. Rev. A 5, 

1852 (1971). [15] 

[4] L. D. Pustilnikov, Theor. Math. Phys. 57, 1035 (1983); [16] 

L. D. Pustilnikov, Sov. Math. Dokl. 35(1), 88 (1987); L. [17] 

D. Pustilnikov, Russ. Acad. Sci. Sb. Math. 82(1), 231 

(1995). 
[5] A. J. Lichtenberg, M.A. Lieberman and R. H. Cohen [18] 

Physica D 1, 291 (1980). 
[6] G. Karner, J. Stat. Phys. 77, 867 (1994). [19] 

[7] S. T. Dembinski, A. J. Makowski and P. Peplowski, Phys. [20] 

Rev. Lett. 70, 1093 (1993). [21] 

[8] J. V. Jose and R. Cordery, Phys. Rev. Lett. 56, 290 

(1986). [22] 

[9] Z. J. Kowalik, M. Franaszek and P. Pieranski, Phys. Rev. 

A 37, 4016 (1988). [23] 

[10] S. Warr, W. Cooke, R. C. BaU and J. M. Huntley, Physica 

A 231, 551 (1996). 
[11] N. Saito, H. Hirooka, J. Ford, F. Vivaldi and G. H. 

Walker, Physica D 5, 273 (1982). 
[12] E. Canale, R. Markarian, S. O. Kamphorst and S. P. de 

Carvalho, Physica D 115, 189 (1998). 



A. Loskutov, A. B. Ryabov and L. G. Akinshin, J. Phys. 
A: Math. Gen. 33, 7973 (2000). 

A. -L. Barabasi, H. E. Stanley, Fractal Concepts in Sur- 
face Crowth (Cambridge University Press, Cambridge, 
1985). 

J. L. Mateos and J. V. Jose, Physica A 257, 434 (1998). 
J. L. Mateos, Phys. Lett. A 256, 113 (1999). 
G. A. Luna-Acosta, G. Orellana-Rivadeneyra, A. 
Mendoza-Galvan and C. Jung, Chaos, Solitons and Frac- 
tals 12, 349 (2001). 

E. D. Leonel and J. K. L. da Silva, Physica A 323, 181 
(2003). 

M. V. Berry, Eur. J. Phys. 2, 91 (1981). 
M. Robnik and M. V. Berry, J. Phys. A 18, 1361 (1985). 
The simplified FUM was introduced by Lichtenberg and 
Lieberman and can be found e.g. in Ref. |3|. 
E. D. Leonel, J. K. L. da Silva and S. O. Kamphorst, 
Physica A 331, 435 (2004). 

We have used the same notation as the original work 0] 
although, at that time, stochasticity was frequently re- 
ferred as to chaotic behaviour. However we stress that the 
transition is from locally to globally chaotic behaviour. In 
the FUM, we mean by locally that the chaotic behaviour 
is confined by two different invariant spanning curves. 



